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Abstract. There exists two types of semi-direct products between a Lie group 
G and a vector space V . The left semi-direct product, G x V, can be con- 
structed when G is equipped with a left action on V. Similarly, the right 
semi-direct product, G N V, can be constructed when G is equipped with a 
right action on V. In this paper, we will construct a new type of semi-direct 
product, G X V . which can be seen as the 'sum' of right and left semi-direct 
products. We then proceed to the parallel existing semi-direct product Eulcr- 
Poincare theory. We find that the group multiplication, the Lie bracket, and 
the diamond operator can each be seen as a sum of the associated concepts 
in right and left semi-direct product theory. Finally, we conclude with a toy 
example and the group of 2-jets of diffeomorphisms above a fixed point. This 
final example has potential use in the creation of particle methods for problems 
on diffeomorphism groups. 

1. Introduction 

Let G be a Lie group and V be a vector space on which G acts by a left action. 
Given these ingredients, we may form the Lie group G k V, which is isomorphic to 
G x V as a set, but equipped with the composition 

(g,v) -x (h,w) = (g ■ h,g ■ w + v) , V(p, v), (h, w) G G x V. 

A standard example of a system which evolves on a left semi-direct product is the 
heavy top, where G = S0(3) and V = M 3 . In contrast, if G acts on V by a right 
action, we may form the right semi-direct product GxV defined by the composition 

(g,v) -x, (h, w) = (g ■ h, w + v ■ h). 

A standard example of such a system is a fluid with a vector-valued advected 
parameter [7]. In any case, it seems natural to surmise that the composition law 

(1) (g,v) -ixi (h,w) = {g-h,g-w + v-h) 

yields a new type of semi-direct product. The first result of this article is that 
([T]) is a valid composition law in some circumstances, where the resulting group is 
dubbed a centered semi- direct product. 

The second result is that the 2nd order Taylor expansions (or 2-jets) of diffeomor- 
phisms over a fixed point is a centered semi-direct product. The main motivation 
behind understanding this example is to allow us to develop particle-based methods 
for fluid simulation, with applications ranging from medical imaging to simulation 
of complex fluids. 

1.1. Background. The semi-direct product is a standard tool used in the con- 
struction of new Lie groups and plays an interesting role in geometric mechanics 
when the normal subgroup is interpreted as an advected parameter. A standard 
example is the modeling of the 'heavy-top', wherein the configuration space is the 
left semi-direct product SO(3) K K 3 [7]. Another standard example is the modeling 
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of liquid crystals, in which we consider the right semi-direct product SDiff(M) x V. 
In this case, SDiff(M) is the set of volume-preserving diffeomorphisms of a volume 
manifold M, and V is a vector space of Lie algebra valued one-forms on M upon 
which SDiff(M) acts by pullback [51 S]. Of course, the tangent bundle of a Lie 
group, TG, is isomorphic to a left semi-direct product G k fj by left-trivializing 
the group structure of TG. Additionally, TG is isomorphic to a right semi-direct 
product G xi g when the group structure of TG is right trivialized [TJ see 55.3]. 
Thus, we see that this method of constructing groups can be found in a number 
of instances. In this article, we introduce a new type of semi-direct product which 
extends the existing semi-direct product theory. 

A motivating example will be a desire to understand the Jet-groupoid of a man- 
ifold M. In particular, it is known that the jet-functor sends the diffcomorphism 
group to a groupoid known as the Jet groupoid O see 512]. As will be illustrated 
in 5jl.4| the isotropy groups of the jet groupoid for 2-jets have a group structure 
which can be written as a centered semi-direct product. A thorough understanding 
of the Jet groupoid can be useful for the creation of new particle-based methods 
wherein the particles carry jet data in addition to position and velocity data. The 
advantage of such a particle method is the possibility for a discrete form of Kelvin's 
circulation theorem [8j. Building such particle methods can be useful in scenarios 
in which one desires to work with the material representation of a fluid. This oc- 
curs in the image registration method known as 'Large Deformation Diffcomorphic 
Matching,' which is used in the field of computational anatomy |3j. Moreover, the 
potential energy associated with the advected parameters in complex fluids is often 
a function of certain gradients which require jet data in order to be advected by a 
diffeomorphism. Thus, keeping track of jet data may play a significant role in the 
construction of particle-based variational integrators. 

1.2. Main Contributions. In this paper, we accomplish a sequence of goals, each 
building upon the previous. In particular: 

(1) In section [2j we define a new type of semi-direct product that we dub a 
centered semi-direct product. 

(2) In proposition |2.2| we derive the Lie algebra of a centered semi-direct prod- 
uct and its associated structures. 

(3) In section [3j we develop the Euler-Poincare theory of centered semi-direct 
product in parallel with the existing theory of semi-direct product reduction 

m- 

(4) In section |4j we describe the centered semi-direct product Euler-Poincare 
equations for a few examples. We present one toy example before presenting 
the theory for an isotropy group of the 2- Jet groupoid. 

Combined, these items allow for a computationally tractable algebraic under- 
standing of 2- Jets and perhaps open the door to applications which were previously 
overlooked by geometric mechanicians. 
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Severo Ochoa Project SEV-2011-0087 and IRSES-project "Geomech-246981". L.C 
owes additional thanks to CSIC and the JAE program for a JAE-Pre grant. The 
work of H.O.J, was supported by European Research Council Advanced Grant 
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1.4. A motivating example. Let Diff(M) denote the diffeomorphisms group of 
a manifold M. For a fixed x € M we may define the isotropy subgroup 

Iso(x) = {ip e Diff(M) | <p(x) = x}. 

Let tp € Iso(a;) and note that T x ip is a linear automorphism of the vector-space 
T X M. In particular: 

Proposition 1.1. The functor "T x " is a group homomorphism from Iso(a;) to 
GL(T X M). 

Proof. Clearly Iso(x) and GL(T X M) are both Lie groups. Let ip,ip € Iso(x). Then 

T x <poT x Tp = T x (<poip). ' ' □ 

This observation has implications for computation for the following reason: By 
definition, T x ip approximates ip in a neighborhood of x € M. Thus, if one desired 
to model a continuum with activity at x, then T x ip carries some of the crucial data 
to do this task. In particular, this is computationally tractable as the dimension of 
GL(T X M) is equal to dim(M) 2 . If dim(M) = n then GL(T X M) = GL(n) is a Lie 
group of dimension n 2 . 



@g£GL(2) 





Figure 1. Depicted is a diffeomorphism with a trivial 2-jet (i.e. 
a linear transformation) and diffeomorphism with a nontrivial 2- 
jet. 

However, the group GL(n) only captures the linearization of a diffeomorphism. 
If we desire to capture some of the nonlinearity then we might consider looking into 
the second jet of these diffeomorphisms (see figure [T]). We can do so by considering 
the functor TT X . Let ip £ lso(x) so that TT x tp is a map from T{T X M) to T(T X M). 
However, T X M is a vector-space so that T(T X M) T X M x T X M . The second 
component represents the vertical component and the isomorphism between TT X M 
and T X M x T X M is given by the vertical lift 



v t (vi,v 2 ) = ^ 



(vi +ev 2 )- 

e=0 



We can therefore represent TT x ip as (T x <p, A v ) where A v : T X M x T X M — > T X M is 
the symmetric (1,2) tensor 

where ip k is the fcth component of ip. In other words, we have the 1-1 correspondence 

TT x p^{A ll A 2 ) 
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where Ax — T x ip and A 2 is given by (|2|). If we denote the set of symmetric T X M- 
valued 2-tensors on T X M by S 2 (x), then this correspondence is given by a map 

* : J 2 \ x x (Diff(M)) -> GL(T.A/) x ^(z) 

where J 2 |* (Diff (M)) is the group of second order taylor expansions about x of 
diffeomorphisms which send x to itself (these are called 2-jets). This allows us to 
write the Lie group structure of J 2 \* x (Diff(M)) as a type of semi-direct product. 
In particular: 



Proposition 1.2. If we represent TT x ip andTT x ip as (Ax,A 2 ) o,nd (Bx,B 2 ) where 

Al =T xl p,B 2 =T x ^A 2 = 0^, 
ip) is given by the composition 



andB 2 = g x ig xJ , then TT x tpo TT x ip = TT x (tpo 



{A 1 ,A 2 ) o (B U B 2 ) = {A x oB 1 ,A 1 oB 2 +A 2 o (Bx x B{)). 
Proof. We find that 



dx, 



dip dip 
dxi dxj 



o ip 



and the second derivative is 
d d 



dxj dx. 



d ( d(p k dip 1 
dxj \ dxi dxi 
d 2 ip k dip 1 dip 



o ip 



dip k d 2 ip l 



dxidx m dxi dx-j dxi dxidx. 



o ip 



A x = 
B x 



Noting that ip(x) — x we can set 

dip k 
dxi 
dip k 
dxi 

and rewrite the equations in the form 
d_ 

d d 



A 2 = 



B 2 = 



d 2 ip k 



dxidxj 
d 2 ip k 



dxidx 1 



dxj dx, 



■(ip k oi/j) 
■((p k oifj) 



Ax-Bx 
Ax-B 2 



A 2 o{Bx x Bx) 



Therefore if we define the composition 

(Ax,A 2 )-(Bx,B 2 ) := ( Ax ■ Bx , Ax ■ B 2 



A 2 o{Bx x Bx)) 



on the manifold GL(T x M)x\/ 2 (T X M; T X M), then * : J 2 \ 
S 2 is a Lie group isomorphism by construction. 



(Diff(M)) 



GL{T x M)x 
□ 



We see that the composition law of Proposition |1.1| is of the form described in 
equation ([TJ. In this paper, we will condense the composition law for 2-jets to 
the algebraic level and study ([I]) in the abstract Lie group setting. Of course, one 
would naturally like to consider diffeomorphisms which are not contained in Iso(x). 
However, this extension brings us into the realm of Lie groupoid theory and will 
need to be addressed in future work. 



2. A CENTERED SEMI-DIRECT PRODUCT THEORY 

In this section, we will discover a new type of semi-direct product. We will 
outline the necessary ingredients for the construction of such a Lie group and we 
will derive the corresponding structures on the Lie algebra. 
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2.1. Preliminary material on Lie groups. Let G be a Lie group with identity 
e € G and Lie algebra g. In this subsection we will establish notation and recall 
relevant notions related to Lie groups and Lie algebras. 

2.1.1. Group actions: Let V be a vector space. A left action of G on V is a smooth 
map pi, : G X V — > V for which: 

Pl(b,v) = v and p L (g,p L (h, «)) = p L (gh,v) , V#, /i € G, Vu e V. 

As using the symbol '/?/_,' can become cumbersome and since we will only need 
a one left Lie group action in a given context, we will opt to use the notation 
g ■ v :— PL(g,v). Finally, the induced infinitesimal left action of g on V is 



tfe 



e=0 



exp(e • £) • v , V£ € g, u e V. 



Similarly, a right action of G on V is the smooth map pn ■ : V x G — > for which: 

pji(«,e) = u and p L (p L (v,g),h) = p L (v,gh) , V#, h £ G, Vu e V. 

Again, we will primarily use the notation v ■ g := pn(v,g) for right actions. The 
induced infinitesimal right action of g on V is given by 



d_ 
~dc 



e=0 



v ■ exp(e • ^) , € g, u e V 



Lastly, we say that the left action and the right action commute if 

(g-v)-h = g-(v-h) 

for any g,h g G and weV. 

2.1.2. Adjoint and coadjoint operators: In this section we will recall the "AD, Ad, ad"- 
notation used in [6] . For jeGwe define the inner automorphism AD : G x G — > G 
as AD(<7, ft) = AD 9 (/i) = ghg^ 1 . Differentiating AD with respect to the second 
argument along curves through the identity produces the Adjoint representation of 
G on g denoted Ad :Gxg^g and given by 



Ad 9 (77) 



d_ 
de 



(AD g (exp(er ? ))) = g ■ 77 • g \ 



for g E G and £ G g. Differentiating Ad with respect to the first argument along 
curves through the identity produces the adjoint operator ad : g x g — > g given by 



ad s (ry) 



d 
1c 



(Ad exp(e?) (r/)) =£.j 7 -t 7 .£. 



£=o 



The ad-map is an alternative notation for the Lie bracket of g in the sense that 

ad(£,7y) ee ad 5 (?7) ee [£,n]. 

For each £ £ g the map ad^ : g — > g is linear and therefore has a formal dua 
ad^ : g* — > g* which we call the coadjoint operator. Explicitly, adj is defined by the 
relation 



(3) 



(ad|(/i),?7) = (p,,ad 6 (rj)) 



for each 77 G g and //eg*. 
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2.2. Centered semi-direct products. In this subsection, we will construct a 
semi-direct product which can be thought of as a 'sum' of a right semi-direct product 
and a left semi-direct product. 

Proposition 2.1. Let G be a Lie group which acts on a vector-space V via left and 
right group actions. Then, the product G x V with the composition law 

(4) (9uVi) ■ (g2,v 2 ) ■■= (fll02,0i ■v 2 + v 1 ■ g 2 ) 

is a Lie group if and only if the left and right actions of G commute. 

Proof. It is clear that G x V is a smooth manifold and that the composition law [4] 
is a smooth map. We must prove that this composition makes G x V a group. 

• That the composition map Q produces another element of G x V can be 
observed directly. Thus 'closure' is satisfied. 

• The identity element is given by (e, 0) € G x V where e € G is the identity 
of G. 

• The inverse element of an arbitrary (g, v) £ GxV is {g , —g vg ) where 
<7 _1 is the inverse of g € G. 

• Given three elements of G x V we find 

(01, «l) ' ((02,^2) • (03, V3)) = (01, «l) • (0203,02 • v 3 + v 2 ■ g 3 ) 

= (010203,01 ' (02 • v 3 + v 2 ■ g 3 ) + vi ■ (.9203)) 

= ((0102)03, (0102) ■ V 3 +01 • («2 • 03) + ("I ' 02) ' 03) ■ 

By the commutativity of the group actions we may equate the above line 
with: 

= ((0102)03, (0102) ■ v 3 + (gt ■ 1)2 ) • 03 + (vi ■ 02) ■ 03) 
= ((0102)03, (0102) • v 3 + (#i ■ v 2 +vi ■ g 2 ) ■ g 3 ) 
= ((0102), 0i • v 2 +vi ■ g 2 ) ■ (03,^3) 

= ((01, «l) • (02,«2)) • (03, «3)- 

Thus, the associative property is satisfied. 

Moreover, all maps in sight including the inverse map are smooth. In conclusion we 
see that GxV with the composition [4] defines a Lie group. Moreover, if the left and 
right actions of G on V do not commute, then we can observe that associativity is 
violated. □ 

Definition 2.1. Given commuting left and right representations of a group G on 
a vector space V , the Lie group GxV with the composition Q is denoted G M V 
and called the centered semi-direct product of G and V. 

It customary to denote the left semi-direct product using the symbol k and the 
right semi-direct product via the symbol x . We justify our use of the symbol M 
in that the concept of centered semi-direct product is merely a 'sum' of a left and 
a right semi-direct product. The formula 1X1= x + x can be used as a heuristic 
throughout the paper. In particular, this heuristic applies to the Lie algebra. 

Proposition 2.2. Let G X V be a centered-semi direct product Lie group. The Lie 
algebra q X V is given by the set q x V with the Lie bracket 

( 5 ) [(£l, u l), (6)V2)]n = (Kl,6] S , (6 ■ «2 + Vl ■ £ 2 ) - (6 -V1+V2- £l)) , 

for £1,62 G 0, V U V2 G V. 
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Proof. Firstly, it is simple to verify that the tangent space at the identity, (e, 0) G 
G x V, is g x V. To derive the Lie bracket, we will derive the the ad-map via the 
Ad and AD-maps. For (g, v),(h,w) G G X V we find 



AD 



(a-h) 



(h,w) = (gh,v ■ h + g ■ w) ■ (g 1 ,-g x -v-g x ) 

= (AB g (h),v ■ hg- 1 +g-wg- 1 - AD g (h) ■ v ■ g' 1 ). 



If we substitute (h,w) with the e-dependent curve (exp(e • ^ 2 ),£ ■ vi) we can 
calculate the adjoint operator, Ad : (G IX V) x (q X V) Q X V. Given by 



Ad (ff ^)(6,u 2 ) 



d 
de 



AD (Sj „)(exp(e ■£ 1 ),e-vx) 



e=0 



(Ad fl (e 2 );»-6jfl 1 +.9-w 2 -.9 1 - Ad 9 (6) -v-g x ). 



If we substitute (g, v) with the i-dependent curve (exp(i£i), ii^) we can differentiate 
with respect to t to produce the adjoint operator ad : (g X V") x (g 1X1 V) — >• g X V. 
Specifically, the adjoint operator is given by 



ad(£i,T>i)(62 ? «2) 



) (Ad (oX p(i. Cl ) it . t , 1 )(^2,W2)) 

(5C2.9 _1 ,f • 6.9 _1 - ff^ff -1 ■ v ■ g- 1 + g ■ v 2 ■ g^ 1 ) 



d 

dt 
d 

~ dt 

= ( ad ?i (6), £l • «2 + «1 • 6 - 6 • «1 - «2 ■ fl) 

= ([Ci,6] B , (6 • «2 + «i ■ 6) - (6 ■ «i + w 2 ■ &)). 

Noting that the ad-map is merely an alternative notation for the Lie bracket 
completes the proof. □ 

We complete this section by defining operations designed to express interaction 
terms between momenta in V in and momenta in G in mechanical systems. 

Definition 2.2. The heart operator Z> : fj x V* — >• V* is defined by 

(6) (£<?a,w)y := (a,£ ■ v - v ■ £) v . 
The diamond operator, <C> : V x V* — > q* , is defined as 

(7) (u<>a,f) fl := (a,v-£-£ -v) B . 

The diamond operator can be seen as the sum of a diamond operator of a left 
semi-direct product and that of a right semi-direct product [7\ . The heart operator 
will allow us to express how momenta in V impact motion in G, while the diamond 
operator will allow us to express the converse. 



3. EULER-POINCARE THEORY 

The Euler-Lagrange equations on a Lie group, G can be expressed by a vector 
field over TG. If the Lagrangian is G-invariant then the equations of motion are G- 
invariant as well and the evolution equations can be reduced. While the unreduced 
system evolves by the Euler-Lagrange equations on TG, the reduced dynamics 
evolve on the quotient TG/G. However, TG/G is just an alternative description 
of the Lie algebra § and so the reduced equations of motion can be described on 
§ where we call them the Euler-P oincare equations. This reduction procedure is 
summarized by the commutative diagram: 
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flow by 'EL' 
TG > TG 



/G 



/G 



flow by 'EP' 

To be even more specific. A Lagrangian L : TG — > E is said to be (right) 
G-invariant if 

L((g,9) ■ h) = L(g,g) 
for all h € G. If £ is G-invariant, then L is uniquely specified by its restriction 
I = L\~ : g — >• M. The Euler-Poincare theorem states the Euler Lagrange equations 

d fdL\ dL (| 



rfi V <9<^ / dg 

on TG are equivalent the Euler-Poincare equations and reconstruction formula 

d (di\ (dl\ - , __ x 

A review of Euler-Poincare reduction is given in (TUJ Ch 13] while a specialization 
to the case of semidirect products with advected parameters is described in [7] . In 
this section we will specialize the Euler-Poincare theorem to the case of centered 
semi-direct products by setting G = G M V. 

To begin let us compute how variations of curves in the group induce variations 
on the trivializations of the velocities to the Lie algebra. Studying such varia- 
tions will allow us to transfer the variational principles on the group to variational 
principles on the Lie algebra. 

Proposition 3.1. Let G M V be a centered semi-direct product and consider a 
curve (g,v)(t) eGNK. Let (£,(<),&(*)) := ($(*)>*(*)) ' G ^ V be 

the right trivialization of (g,v)(t). An arbitrary variation of (g,v)(t) is given by 

(Sg,Sv)(t) = (v g ,Vv)(t) ■ (g,v)(t) G T {g . v){t) {G\AV), 

where {r)g,r) v )(t) £ g 1X1 V . Given such a variation, the induced variation on (£ g ,£ v ) 
is given by 

(8) (S€g,S€v) = (Vg - ad ? 9 Vg,Vv + (Vg€v + Vvtg) ~ (^gVv + ZvVg)) 

d 



Proof. For any Lie group, G, and any curve g(i) € G, the variation of £(t) := 
g(t) induced by the variation 6g(t) — fj(t) -git) is 5£ — fj— [£, fj]. For matrix 

groups see (TU1 Theorem 13.5.3] and [2] for the general case. If we set G = G ixi V 
and use the bracket derived in Proposition |2.2| then the theorem follows. □ 



Now that we understand the relationship between variations of curves in G IX V 
and the induced variations in g IX V we can state the Euler-Poincare theorem for 
centered semi-direct products. 

Theorem 3.1. Let L : G IX V ->• K be (right) G IX V -invariant, and let I : g IX V -> 

K 6e ite reduced Lagrangian. Let (g,v)(t) e G IX and denote the right trivialized 
velocity by (£g f := (g,v)(t) ■ (g, v)^)^ 1 . Then the following statements are 

equivalent: 
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(i) Hamilton's principle holds. That is, 
(9) <5 / ' L(g(t),g(t),v(t))dt = 



for variations of (g,v)(t) with fixed endpoints. 
(ii) (g,v)(t) satisfies the Euler- Lagrange equations for L. 
(Hi) The constrained variational principle 

(10) s f 1 i(Z a (t),z v (t))dt = o 

J to 

holds on g x V for variations of the form 

(11) (5£ g , 5£ v ) = {rig - ad Cs rj g , »)„ + r) g £ v - £ v r) g + r, v £, g - £, g r, v ). 

where {ri g ,ri v )(t) is an arbitrary curve in g 1X1 V which vanishes at the end- 
points. 

(iv) The Euler- Poincare equations 

d ( 5i\ ,« / 5£ \ , A 51 



d 

It \<K, I V " <Kr 

hold on gMV . 



Proof. The equivalence (i) and (ii) holds for any configuration manifold and so, 
in particular it holds in this case. 

Next we show the equivalence (Hi) and (iv). We compute the variations of the 
action integral to be 

= it ^W^ 9 ~ ^ V9 ) + (sY ,f>V+T,3 ^' V ~ ^a+Vvtg - ZgVv) 
and applying integration by parts and equation (|3| we find 



51 

dt \IK 



— adt 



5t_ 

Wg 



d dl 

dt d£ v 



i 81 
/ dl 

+ \Wg^ 



\ I Si 

VgZv ~ CvVg ) + {jT'VvZg ~ ZgVv fdt 



5$ v 



+ 



\5£ v 



51 

dt \5tfg 

8^ 
dt \5£ v 



81 

Wg 



5£ v 



, Vg 



By noting that (ji g , ri v )(t) is arbitrary on the interior of the integration domain, the 
result follows. 

Finally, we show that (i) and (Hi) are equivalent. The G— invariance of L implies 
that the integrands in (|9| and (10) are equal. However, by Proposition 3.1 all the 
variations of (g, v) (t) with fixed endpoints induce, and are induced by, variations 
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(8£g,8£v)(t) <E 1X1 V of the form given in equation (11). Conversely if (i) holds 
with respect to arbitrary variations (8g,8v), we define 



(Vg,Vv)(t) = (8g,8v) ■ (g,v)' 



to produce the variation of £„) given in equation ( 11 ). 



□ 



Remark 3.1. There is a left invariant version of theorem (3.1 1 in which (£ g ,£ v ) : = 
• (<7,t>) and L is left G 1X1 V -invariant. In this case the Euler-Poincare 
equations take the form 



dt \dt[g 



ad 



d f 81 
dt \d^, 



81 
o£,v 



S£^ 
8jv 



= 0. 



Remark 3.2. There is a version of semi-direct product mechanics wherein the 
vector-space V is a set of advected parameters as in . In this case we impose 
the holonomic constraint 

v = g ■ v + v ■ g 
and the set of admissible variations in g X V become 

8£g =Vg- [Zg, Vg] , 8v = rj g ■ v + v ■ T) g . 

If we do this, the ty-term is removed and |^ equation is replaced with a holonomic 
constraint. In particular we find that 



d f Si 



± ad^ 



+ a 81 



dv 



8£g 

= £g-V + V-Z g . 



where we use a plus sign for right trivialization and a minus sign for left trivializa- 
tion. 



4. Examples 

In this section we will present two examples of Euler-Poincare equations on 
centered semidirect products. This first is a toy example designed to illustrate how 
computations of the diamond and heart operators can be done in practice. The 
second example is designed to 2- Jets as described in subsection [E4] 

4.1. A toy example. Consider the group GL(n) and let Mat(n) denote the vector 
space ofnxn real matrices. Noting that GL(n) acts on Mat(n) by left and right 
multiplication, we can define the composition law on the Lie group GL(n) X Mat(n) 
by: 

(A, v) ■ (B, w) = (AB, Aw + vB). 

Moreover, we can identify fl(*(n) with gl(n) and Mat(n)* with Mat(n) by the matrix 
trace pairing (A,B) — Tr(A T B). This allows us to calculate the heart operator 
9 : fll(n) x Mat(n)* -> Mat(n) as 
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(AVw,v) = {w, A ■ v - v ■ A) 

= trace (w T (A ■ v — v ■ A)) 
= trace (w T ■ (A ■ v) — w T (v ■ A)) 
— trace ((w T ■ A)v — (A ■ w T ) ■ v) 
= trace ((w T ■ A — A ■ w T ) ■ v) 
= trace [(A T w - w ■ A T ) T ■ v) 
= (A T w - wA T , v) 

Therefore, 

AVw = A T w - wA T . 
By a similar calculation, diamond operator is found to be 

v()w = v T w — wv T , 
and the coadjoint action on GL(n) is given by 

&d* A (aA) — A T ■ ola — oiA • A T . 

Now, we have all the ingredients to write the Euler-Poincare equations. Given a 
reduced Lagrangian I : gl(n) X Mat(n) ->lwe may denote the reduced momenta 
by 

St Si 

^ = s~e 1 = Tv- 

where (£, i>) € Qi(n) M Mat(n). The Euler-Poincare equations can be written as 
fi = - fi£, T ) + w T 7 - 7« T 

7 = C T 7~7C T - 

4.2. An isotropy group of a 2- Jet groupoid. In proposition 1 1 . 1 1 we illustrated 
how the set of 2-jets of diffeomorphisms of the stabilizer group of a point x E M 
is identifiable with a centered semidirect product. In particular, if dim(M) = n we 
can consider the group GL(n) M 5^, where 5^ is the set of (1, 2)-tensors which are 
symmetric in the covariant part. For the moment we shall consider the larger space 
of all (1, 2)-tensors denoted T^ 1 . If we let ei, . . . , e„ G M. n be a basis with dual basis 
e 1 , . . . , e™ e (E™)* we can write an arbitrary element of 7j as 

T = T} k ei®e j ®e k . 

The left action of GL(n) on T^ 1 is 

g-T:= T] k {g ■ e,) ® ® e fe = T]J^ X ® e j ® e fc 

while the right action is 

T ■ g := 7™ fc e 8 ® ( 5 T • e*) ® ( 5 T • e fc ). 

Clearly these actions commute, and so we may form the centered semidirect product 
Lie group GL(n) X T2 1 . 

Let us now focus on the Lie algebra. The Lie algebra gi(n) is equivalent to 
and the Lie bracket is then given in the bases ej ® e- 7 by 

where £ = £je, ® e J and r\ — ff^i ® e J . We can use the dual basis e 1 ® ej to sec 
that the coadjoint action of £ on /j, = /.ije 1 ® e,,- is given by 

ad^M=(Mfc^-M^)e < ®e j . 



12 



LEONARDO COLOMBO & HENRY O. JACOBS 



By differentiation we see that the infinitesimal left and right actions of Ql(n) on 7a 
are given by 

£ • T = Tj k ^ei ® ® e fc 



T • £ = 7] 



e 4 <g> (£f • e z ) <g> e k + e t <g> e 1 <E> (ff ■ e l ) 
= (T i Uj+^e|)e i ®^®e fe . 

If we choose an arbitrary element a € (T2)* = TJ 2 given by 



a = e* ® e,- ® e fc 



we find that 

(a,t-T) = (al k &TL = (a?T? k )Q 



•pi / 

Therefore the heart operator is given by 

= (£af - - aft,V ® ej ® e fc 

and the diamond operator is 

o<>T = (a?*2& + - af 2* )e< ® ej . 

Given a reduced Lagrangian £ : fll(n) M 7^ — > K we can denote (J, = J| and 7 = 
In terms of the basis e l ®e,,- and ej®e ? '®e fe we may write the (right) Euler-Poincare 
equations as: 



By restricting T^ 1 to the subspace S\, we can obtain a Lie group which models 2-jets 
of diffeomorphisms as demonstrated in proposition 1.2 This example provides a 
first step towards the creation of higher-order, spatially accurate particle methods 
[H section 4.2]. Moreover, the data of 2-jets is necessary for the advection of 
quantities seen in complex fluids in which the advected parameters depend on 
gradients of the flow [H [5] . Therefore, the structures described here may prove 
useful in the construction of particle-based integrators for complex fluids as well. 



5. Conclusion 

In this paper, we have presented a variant of traditional semi-direct products, 
dubbed centered semi-direct products, and we have illustrated the associated Euler- 
Poincare theory. The diamond operator, the group multiplication, and the Lie 
bracket can all be seen as sums of the associated concepts for left and right semi- 
direct products. As a result, the Euler-Poincare theory associated with centered 
semi-direct products can also be seen as a sum of the left and right invariant Euler- 
Poincare theories for semi-direct products. Presently, many of these constructions 
remain fairly theoretical. However, an isotropy group of the groupoid of 2-jets of 
diffeomorphisms of a manifold can be seen as a centered semi-direct product. This 
has potential applications in simulation of complex fluids. We hope this paper 
provides a stepping stone towards realizing this application. 
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